Correlation tuned cross-over between thermal and nonthermal states following 

ultrafast transient pumping 
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We examine electron-electron mediated relaxation following excitation of a correlated system 
by an ultrafast electric field pump pulse. The results reveal a dichotomy in the temporal evolu- 
tion as one tunes through a Mott metal-to-insulator transition: in the metallic regime relaxation 
can be characterized by evolution toward a steady-state electronic distribution well described by 
Fermi-Dirac statistics with an increased effective temperature; however, in the insulating regime 
this quasithermal paradigm breaks down with relaxation toward a nonthermal state with a more 
complicated electronic distribution that does not vary monotonically as a function of energy. We 
characterize the behavior by studying changes in the energy, photoemission response, and electronic 
distribution as functions of time. Qualitatively these results should be observable on short enough 
time scales that the electrons behave like an isolated system not in contact with additional degrees 
of freedom which can act as a thermal bath. Importantly, proper modeling used to analyze ex- 
perimental findings should account for this behavior, especially when using strong driving fields or 
studying materials whose physics may manifest the effects of strong correlations. 

PACS numbers: 78.47. J-, 71. 10.Fd,79.60.-i 



Optical reflectivity, pj-foj photoemission spectroscopy, 
[E ErEa] and resonant x-ray scattering 16- 18[ are equi- 
librium methods which in the time domain are ideally 
suited to studying dynamics of novel ordered phases or 
collective excitations. @-H, @-[l(J EMI 0n sufficiently 
short time scales, the initial recovery in these systems 
following an ultrafast pump pulse should be dominated 
by electron-electron scattering which on its own can drive 
the system into a new steady-state. Conventional anal- 
ysis has been based on a quasithermal paradigm ( "hot 
electron" or multi-temperature models) ;[l9l.[20j however, 
there have been few tests of the validity of its underlying 
assumptions as a function of the strength of electronic 
correlations, [2l| in particular as one tunes between the 
two regimes of a metal-to-insulator transition (MIT).[22[ 

The MIT driven by electronic correlations usually is ac- 
companied by a number of interesting ordering phenom- 
ena among the spin, charge, and orbital degrees of free- 
dom in a material. An understanding of the key physics 
which leads to these emergent phases is often at the 
heart of pump-probe experiments in condensed matter 
systems, including high-T c cuprate superconductors, [23 



what can be learned by tuning the interaction parameters 
of model systems simulated in fermionic or bosonic cold 
atom mixtures and performing the experimental equiva- 
lent of time-resolved, pump-probe measurements. 27- 2fjj| 
In this Letter we discuss the evolution of an electronic 
system described in equilibrium by a simple Hamiltonian 
representing a correlated electronic system possessing a 
MIT tuned by the strength of the Coulomb repulsion U. 
In order to avoid approximate treatments of either cor- 
relations or applied fields, we chose to study the spinless 
Falicov-Kimball model[30( whose effective Hamiltonian is 
given by 
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nickelates, manganites, ruthenates, vanadates. [22] and 
even organic materials. 24-2(| A number of experimen- 
tal parameters can be used to tune across the MIT in- 
cluding doping and chemical substitution, pressure, and 
applied fields. What can be learned about the underly- 
ing physics leading to these phases as a function of these 
key parameters requires an understanding of the proper 
paradigm in which to ask the relevant questions and con- 
duct analysis of experimental data. This is in addition to 



This model describes itinerant conduction electrons hop- 
ping between lattice sites with an energy t* and chemical 
potential \i that experience Coulomb repulsion with an- 
other species of localized electrons with an occupation Wi 
on each site. This model can be tuned through a Mott 
MIT at half-filling by adjusting the electron-electron in- 
teraction U with U c = y/2t*. Here we wish to un- 
derstand, in general terms, how electron dynamics are 
affected by both strong fields and strong correlations; 
therefore, our exact treatment for this model is of more 
general interest than using a materials-specific Hamilto- 
nian that may require a number of approximations to 
affect a full solution in the time-domain. 

We model the transient pump pulse as a spatially uni- 
form, harmonic, electric field with a Gaussian envelope 
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FIG. 1: (a) The electric field pump pulse that leads to 
the temporal evolution shown in panels (b) and (c) has a 
maximum intensity _E max =24i5 (normalized in the panel), 
modulation frequency cj p = 0.5t*, and characteristic width 
<7 P = 5/t*. (b) & (c) Time-resolved pump-probe photoemis- 
sion response for (b) metallic (U = 0.5t*) and (c) insulating 
(U — 2t*) systems determined for a probe pulse with charac- 
teristic width Ob = 2/t* . 



[see the temporal profile shown in Fig. HJa)] incorpo- 
rated via the Peierls' substitution|3l| in the Hamiltonian 
gauge. The temporal evolution is simulated using an ex- 
act, nonequilibrium formulation of DMFT. 32H36| Fol- 
lowing standard convention, the energy unit is taken to 
be t* throughout this work and the standard unit for time 
is 1/t*. Conversion to physical units and the effective 
electric field scale, denoted E a , as well as details about 
the evaluation of an effective quasithermal response, can 
be found in the Supplementary Material. 

Figures Gib) and (c) show the characteristic photoe- 
mission response [37- 39] as a function of time delay for 
representative metallic and insulating systems photoex- 
cited by the pump-pulse shown in Fig.[lja). The width of 
the photoemission probe pulse influences both the tem- 
poral resolution and energy resolution of the resulting 
spectrum, chosen here to strike a balance between the 
two. 

For metallic correlations [Fig H^b), U = 0M*] the 
pump pulse narrows and shifts the response toward the 
equilibrium Fermi level {ui = 0t*). Following the pump 
pulse, we observe a rapid relaxation toward a signifi- 
cantly broader spectral distribution characteristic of a 
new steady-state. For stronger correlations on the insu- 
lating side of the Mott MIT [FigQJc), U = 2t*} the pump 
pulse narrows the response below the equilibrium Fermi 
level and transfers spectral weight across the insulating 
Mott gap (centered at the Fermi level). Similar to the be- 
havior observed for weak correlations, we find a broader 
spectral distribution following the transient pump, al- 
though in this case a significant remnant spectral weight 
transfer across the gap persists, even in the long-time 
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FIG. 2: Instantaneous power delivered to (a) metallic (U = 
0.5t*) and (b) insulating (U = 2t*) systems discussed in 
Fig- HJ ( c ) & (d) As a closed system, the change in total 
energy may be determined from the local propagator (solid 
lines) or integrating the instantaneous power (dashed lines). 
The constant offset (dotted lines) represents the initial, equi- 
librium value in each case. The insets highlight times near 
the center of the pump pulse between —2.5/t* and 2.5/t* that 
show rapid variation in each quantity. 



steady-state. This rapid evolution following the decay of 
the pump pulse can be attributed to relaxation mediated 
by electron-electron scattering. 

An increase in the effective temperature would nat- 
urally lead to broader features and a redistribution of 
weight across the insulating Mott gap, as the spectral 
function in equilibrium is temperature independent for 
this model. The highlighted (red) traces in Fig. [T] can 
be used to assess whether the steady-state is representa- 
tive of the system in equilibrium at an elevated, effective 
temperature. However, as we now will show, our results 
highlight a distinct dichotomy in the temporal evolution 
of as one tunes across the MIT. In the metallic regime the 
quasithermal picture remains valid with electron-electron 
mediated relaxation characterized by evolution toward an 
effectively thermal steady-state. However, tuning corre- 
lations across the MIT causes a breakdown in the qua- 
sithermal paradigm in the insulating regime with clearly 
nonthermal relaxation even at long times. 

We first determine the appropriate temperature based 
on the change in total energy of the system. Figure [2] 
shows the power delivered to the system by the pump 
pulse as well as the time-dependent total energy for the 
systems discussed in Fig.[TJ Figures [2jc) and (d) show a 
comparison between the total energy and the change in 
total energy determined by integrating the instantaneous 
power. As one can see, the simulation properly accounts 
for the energy delivered by the pump pulse with small 
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FIG. 3: Quasithermal fits of the steady-state photoemission 
response (characteristic probe width <7b = 2/t*) for (a) metal- 
lic (U = 0.5t*) and (b) insulating (U = 2t*) systems, respec- 
tively, with the effective temperature shown in parenthesis. 
The fits (solid black) are determined by extracting the tem- 
perature from the total energy, while the best-fit quasither- 
mal response (dashed red) has been determined by the least- 
squares method (LSQ). 



deviations near time 0/t* (near the center of the pulse) 
where the power (and instantaneous current) changes 
rapidly. The offset simply reflects the initial energy in 
equilibrium. Converting the total energy at long times 
to an effective temperature (assuming validity of the "hot 
electron" model) , we find an increase of ~ 8 times for the 
weakly correlated, metallic system and ~ 43 times for the 
strongly correlated, insulating system for identical driv- 
ing fields. 

Figure [3] shows the result of quasithermal fits to the 
photoemission response in both [(a), U = 0.5t*] the 
metallic and [(b), U = 2.0t*] insulating regimes. In each 
case the characteristic probe width is fixed as in Fig. [T] 
In the weakly correlated metal the simulated response 
closely matches that derived from the quasithermal fits 
with only small differences which grow upon increasing 
the pump strength and frequency as would befit simple 
expectations (see the Supplementary Material). In this 
case the "best fit" also has been determined from a simple 
least-squares fit (LSQ), but with a similar conclusion. 

In contrast, the steady-state response in the insulating 
regime shows significant deviation from a quasithermal 
best fit, determined using either method, clearly indicat- 
ing a breakdown in the "hot-electron" model. Issues asso- 
ciated with convergence of our self-consistent method re- 
strict our studies to relatively strong driving fields; how- 



ever, we can infer that similar observations should hold 
as the strength of the transient field is reduced, although 
the deviation from the quasithermal paradigm is likely 
to be noticeable only in the tail of the spectral function 
across the Fermi level, or, equivalently, the extracted, 
effective Fermi-Dirac distribution. Observations of non- 
thermal behavior already have been made experimentally 
for comparatively weaker driving fields. 18 1 

We also examine the instantaneous electronic distribu- 
tion providing a snapshot of the response to the applied 
electric field and subsequent relaxation due to electron- 
electron scattering. Consider a simple noninteracting 
band metal. The initial distribution follows usual Fermi- 
Dirac statistics f(£k — with preferential occupation of 
the lower-energy states according to temperature. When 
driven by an electric field, an electron initially at momen- 
tum k shifts to k — eA(i) (the standard Pcicrls' substitu- 
tion). This behavior can be extracted directly from the 
gauge-invariant lesser Green's function for the system. 

Consider the results shown in Fig. Plotted versus 
the band energy and the normalized electron velocity 
v k = Vk • E, the line dividing the occupied and empty 
states rotates at a rate given by — dA(t) jdt-E = E(i) • E, 
simply the magnitude of the electric field as a function 
of time. One may remove this rotation by going-over to 
an instantaneous frame which corresponds to examining 
the distribution functions in a particular gauge where one 
observes a static Fermi-Dirac distribution with no tem- 
poral dynamics. While the underlying physics remains 
unchanged, visualizing the results depends on whether 
one works with the gauged or gauge-invariant functions. 
We choose the latter (gauge-invariant formulation) as il- 
lustrated for a simple half-filled one-dimensional nonin- 
teracting band metal in Figs. |4ja-d) and used to display 
the instantaneous electronic distribution functions for the 
metallic and insulating regimes shown in Figs. Hie) and 
(f), respectively. 

Incorporating the influence of interactions has a num- 
ber of effects on the temporal behavior. First, the initial 
distribution broadens due to the electron-electron inter- 
actions, but remains independent of «k and monotoni- 
cally decreases as a function of £k- The applied field 
drives particles via the Peierls' substitution, but now 
electron-electron scattering randomizes the distribution 
and eventually quenches Bloch oscillations initially ob- 
servable in the distribution function. After the pump 
pulse decays, electron-electron scattering drives the dis- 
tribution toward a more-or-less static, steady-state, pat- 
tern. If the distribution depends on or, as primarily 
observed in these simulations, is no longer monotonically 
decreasing as a function of £k, then the quasithermal 
paradigm cannot hold. 

Figs. IUe) and (f) show the equal-time distribution 
function for a sequence of times in the metallic and in- 
sulating regimes, respectively. In Fig. Hfe) weak corre- 
lations allow the electric field pulse to easily shift the 
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FIG. 4: (a)-(d) A simple cartoon depicting the influence of 
an applied electric field on a noninteracting one-dimensional 
band metal at half-filling, (a) & (c) The band energy (ek, 
red highlight) and normalized band velocity (v^, blue high- 
light) distributions in equilibrium and under the influence of 
an applied electric field, respectively. In one-dimension each 
band energy is associated with two band velocities (right (+) 
and left (-)). In higher dimensions a distribution of velocities 
is associated with each band energy, (b) & (d) The equiv- 
alent distribution shown in a two-dimensional band energy- 
velocity space. The grey shading is a guide to the eye. (e) 
& (f) Equal-time band energy-velocity distribution functions 
for the metallic (U = 0.5t*) and insulating (U = 2t*) regimes 
of Fig. 1(b) and (c), respectively, for various times (units 
of 1/t*). The band energy- velocity space spans a radius of 
3.9t* with the center of each plot at the origin ( Ot* , Ot*) as 
indicated by the band energy-velocity axes in the top pan- 
els. A time-lapse sequence of images for both cases with finer 
resolution can be found in the Supplementary Material. 

distribution function, similar to the expected behavior 
for a simple band metal. The system maintains a well- 
defined "Fermi edge" , as observed in equilibrium, for all 
but those times with the strongest electric field near the 
center of the pulse (0/t*). The edge reforms as the pump 
pulse decays and at the longest simulation times a signif- 
icantly wider edge appears indicative of the higher ef- 
fective temperatures, and validity of the quasithermal 
paradigm, used to describe the observed photoemission 
response in Fig. [3] 

In Fig. |U[f), electronic correlations on the insulating 



side of the Mott MIT produce an equilibrium distribution 
also with a well-defined "Fermi edge" . These correlations 
provide an initial resistance to the influence of the applied 
pump pulse at short times and significantly scramble the 
electron redistribution at long times. Relaxation through 
electron-electron scattering induces a partial reformation 
of the edge at the longest simulation times; however, the 
system retains a significant nonthermal distribution of 
weight characterized primarily by a nonmonotonic de- 
pendence on £k- A sequence of snapshots for both cases 
with finer time resolution and snapshots for the observed 
behavior with an alternative pump pulse can be found in 
the Supplementary Material. 

These results reveal a dichotomy in the evolution of 
transiently excited electrons as one tunes across the Mott 
MIT. The quasithermal picture, which has served to 
underpin much of the experimental analysis, remains 
valid in the metallic regime where relaxation can be 
characterized by evolution toward an effectively ther- 
mal steady-state. Tuning correlations to the insulat- 
ing regime, across the MIT, causes a breakdown in this 
paradigm as one clearly observes relaxation toward a non- 
thermal state. On short time scales where the electronic 
system behaves like an isolated system not in contact 
with additional degrees of freedom, one must take this 
dichotomy into account when performing experimental 
analysis. Additional interactions, not considered here, 
should dominate the much longer time recovery where the 
system must naturally return to its original equilibrium 
through coupling to the crystal lattice (electron-phonon 
coupling) and eventually ballistic and diffusive transport 
of the delivered pump energy to the material's bulk and 
subsequently the environment. 
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Supplementary Material: Correlation tuned 
cross-over between thermal and nonthermal 
states following ultrafast transient pumping 



SI: Scaling 

Scaling of physical quantities in the nonequilibrium 
DMFT formalism begins by considering the proper en- 
ergy scale via the noninteracting kinetic energy. Con- 
sider only nearest-neighbor hopping on the infinite di- 
mensional hypercubic lattice. Then the band energy 
£k = linid->oo — 2iX)i=i cos(fcj), where d is the lattice di- 
mension. The hopping is chosen to scale si as t — t* /2\fd 
and we take t* as the base energy unit for the prob- 
lem. This scaling leads to a Gaussian noninteracting 
density of states and a hnite effective bandwidth. Since 
the direction of the applied field was chosen along the 
hypercubic body diagonal, the Pierels' substitution nat- 
ural produces a second "band energy" given by — 
lim^oo — (t*/\/(tj sin(fcj), which disappears from 

the problem in the limit where the strength of the applied 
field goes to zero. These two band energies have a Gaus- 
sian joint density of states which appears in the Hilbert 
transform used to obtain the local Green's function in 
the iterative DMFT process. 

With the effective energy scale t*, time is measured 
in units of 1/t* and the magnitude of the effective field 
would have units of t*/e. To convert these to some phys- 
ical scale, first we choose t* ~ 0.1 eV. Our fundamental 
unit of time becomes h/t* ~ 6.6 fs; and if we assume 
the characteristic length scale a ~ 3A, the character- 
istic unit for the field E Q ~ t*/ea ~ 33 mV/A. To 
add some additional perspective, the equilibrium tem- 
perature for all cases considered in the main text would 
be T = O.lt* ~ 116 K; and for the pump pulse in the 
main text w p = 0M* ~ 76 THz, a p = 5/t* ~ 33 fs, 
and £ max ~ 800 mV/A (a large value, but not unrealis- 
tic given the field strengths currently employed in some 
pump-probe experiments) . 



S2: Power and Energy 

We evaluate the instantaneous power, kinetic energy 
(KE), potential energy (PE), and total energy (see Fig. 2 
in the main text and Fig. ISlj) . The instantaneous power 
delivered by the pump pulse can be determined in a 
straightforward manner by first evaluating the instan- 
taneous current 



(j(*)> 



-ei 



X] V k-eA(t)Gk (t,t), 



(SI) 



or alternatively, as we do in this case, one can ex- 
press both the band velocity and momentum-dependent 
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FIG. SI: (a) The kinetic energy (KE) as a function of time 
for the two interaction regimes in the main text, (b) The 
corresponding potential energy (PE) as a function of time. 
The insets highlight those times near the center of the pump 
pulse where one sees rapid variation in the KE as a function of 
time. The effect of the pump pulse manifests as an increase in 
KE in both the metallic and insulating regimes; however, the 
PE shows a distinct difference, as one naively expects, with a 
significant PE increase as a function of time in the insulating 
regime due to the applied pump pulse. 



lesser Green's function in terms of the two band energies 
and convert the sum over momentum into a weighted 
integral over {e,e} as one does in evaluating the lo- 
cal Green's function. Then the instantaneous power 
P(t) = (j(t)) E(t), as both the current and field point 
along the hypercubic body-diagonal in our simulations. 
Since we are dealing with a closed system, the integral of 
the instantaneous power as a function of time gives the 
total energy delivered by the applied field. 

We evaluate the KE in much the same fashion as the 
instantaneous current, S3 since one can merely replace the 
velocity Vk with the band energy and drop the electric 
charge (again one equivalently performs an integral over 
{e,e} in the numerical evaluation). The total system 
energy is determined from the time derivative of the local 
lesser Green's function 83 



Total Energy : 



(t,f) 



Of 



(S2) 



analogous to the equilibrium expression. While we could, 
in principle, evaluate the PE from a product of the self- 
energy and Green's function analogous to the equilibrium 
expression, 83 the system is closed and the PE = Total 
Energy - KE. Each of these quantities has been extrapo- 
lated independently to the limit of zero discretization on 
the Keldysh contour. 
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FIG. S2: The temperature T versus the total system energy 
in the Falicov-Kimball model for the two values of interac- 
tion U studied in the main text. Points highlight the energy- 
temperature values corresponding to the pump pulse from 
the main text (Pump #1) and a second (Pump #2) appear- 
ing only in this Supplementary Material with characteristics 
Emax^OEo, uip = l.Ot* , and a p = 5/t* (see Fig. [S5] for a 
comparison) . 



S3: Quasithermal Fitting 

For a given interaction strength, the local density of 
states N(uj) in the Falicov-Kimball model is temperature 
independent. Therefore, in a quasithermal fit to the pho- 
toemission response, the temperature enters only in the 
effective Fermi-Dirac distribution function which gives 
the effective local lesser Green's function —iGf^ Jw) = 
N(u))f(uj,T). One measure for quasithermal effective 
temperature comes from the total energy (see Fig. 2 in 
the main text). Figure [S2l shows the temperature as a 
function of total energy for the two interaction strengths 
discussed in the main text. Here, the total energy for 
a given temperature has been determined from the first 
moment of N(u>)f(u),T). Points have been included in 
Fig. [S2] for two different pump pulses: the one listed in 
the main text (Pump #1) and a second (Pump #2) dis- 
cussed only in the Supplementary Material. Using these 
effective temperatures, the quasithermal photoemission 
response can be determined straightforwardly. 84 The re- 
sults appear as solid lines in Fig. 3 of the main text or 
Fig. [S4]for the two pump pulses, respectively. 

We also have determined the effective quasithermal re- 
sponse for a large range of temperatures. Figure [S3] shows 
some representative results for a few temperatures com- 
pared to the simulated data. The inset in each panel of 
Fig. [S3] shows the sum of the square difference between 
the quasithermal response and the simulation data as a 
function of temperature. The quasithermal response for 
the temperature corresponding to the minimum xtsQ m 
each case determines the least-squares fit appearing in 
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FIG. S3: The equilibrium photoemission response for [a, 
U = 0.57;*] metallic and [b, U = 2.0t*] insulating regimes dis- 
cussed in the main text for several different temperatures T. 
The insets show the least-squares fit to the simulated response 
function at long times used to determined the least-squares 
"best-fit" response appearing in Fig. 3(a) and (b), respec- 
tively. Black dots represent data from the simulation as in 
the main text (see Fig. 3), shown only for every tenth data 
point for clarity. The full data set has been used to determine 
the least squares fit in each case. 

Fig. 3(a) and (b) of the main text, respectively. A sim- 
ilar analysis has been conducted for the second pump 
pulse (P. #2) appearing in Fig. IS4I 

S4: Movies 

As stated in the caption to Fig. 4 in the main text a 
series of movies with fine-grained time resolution have 
been produced for both the (a) metallic and (b) insulat- 
ing regimes under study subject to the two pump pulses. 
Figure [S5] shows the temporal profile of the electric field 
E for the two cases. The notation used for the movies is 
SM followed by a numeral representing the pump number 
and ?7i or i representing the metallic or insulating regime, 
respectively, e.g. the movie following the time evolution 
of the equal-time distribution function for the metallic 
system subject to the pump discussed in the main text 
(Pump #1) is SMlm and that following the evolution for 
the insulating system subject to the second pump (Pump 
#2) is SM2L As in the main text, the band energy- 
velocity space spans a radius of 3.9t* with the center at 
the origin (Or;*, 0i*) as indicated by the axis at the center 
of the plots. The respective time for the frames of each 
movie has been annotated in the upper left hand corner 
in units of 1/t* with the interaction strength U in units 
of t* on the right. In each case, the distribution function 
"rotates" in band energy- velocity space. In the metallic 
regime, the "Fermi edge" present in equilibrium reforms 
at long times with a significantly wider distribution in- 
dicative of the increase in effective temperature; at long 



8 




times one observes the nonmonotonic dependence on 
in the insulating regime indicative of a nonthermal state 
following the ultrafast transient pumping. 

Supplemental References 

51. W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 
324-327 (1989). 

52. J. K. Freericks, Phys. Rev. B 77, 075109 (2008). 

53. V. M. Turkowski and J. K. Freericks, Phys. Rev. B 
73, 075108 (2006). 

54. J. K. Freericks, H. R. Krishnamurthy, and T. Pr- 
uschke, Phys. Rev. Lett. 102, 136401 (2009). 



CO (t* ) CO (t*) 



FIG. S4: Quasithermal fits of the steady-state photoemis- 
sion response (characteristic probe width <Tb = 2/t*) for (a) 
metallic (U = 0.5t*) and (b) insulating (U = 2t*) systems 
subject to Pump #2 (P. #2), respectively, with the effective 
temperature shown in parenthesis. The fits (solid black) are 
determined by extracting the temperature from the total en- 
ergy, while the best-fit quasithermal response (dashed red) 
has determined by the least-squares method as described in 
Section S3. 




FIG. S5: Temporal profile of the electric field pump pulse E 
for two different cases: the one discussed in the main text 
(Pump #1) and the second at higher fundamental frequency 
and field strength (Pump #2) as described in Fig. IS2I 



